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^ ■ SONYC - Substellar Objects in Nearby Young Clusters - is a survey program 

to investigate the frequency and properties of brown dwarfs down to masses 
below the Deuterium burning limit in nearby star forming regions. In this second 
paper, we present results on the ~1 Myr old cluster p Ophiuchi, combining our 
own deep optical and near-infrared imaging using Subaru with photometry from 
the 2-Micron All Sky Survey and the Spitzer Space Telescope. Of the candidates 
selected from iJK s photometry, we have confirmed three - including a new brown 
dwarf with a mass close to the Deuterium limit - as likely cluster members 

OO ■ 

. through low-resolution infrared spectroscopy. We also identify 27 sub-stellar 

candidates with mid-infrared excess consistent with disk emission, of which 16 
are new and 11 are previously spectroscopically confirmed brown dwarfs. The 
high and variable extinction makes it difficult to obtain the complete sub-stellar 
population in this region. However, current data suggest that its ratio of low- 
mass stars to brown dwarfs in similar to those reported for several other clusters, 



though higher than what was found for NGC 1333 in lScholz et al.l (120091 ) . 
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Introduction 



Understanding the origin of the stellar initial mass function (IMF) is one of the major 
topics in astrophysics. The low-mass end of the IMF, in particular, has been subject o f 
numerous observational and theoretical studies over the past decade (see lBonnell et al.ll2007l ). 



SONYC - Substellar Objects in Nearby Young Clusters - is an ongoing project to 
provide a complete census of the brown dwarf (BD) and planetary mass object population 
in nearby young clusters, and to establish the frequency of substellar mass objects as a 
function of cluster environment. The resulting catalog of substellar mass candidates will 
provide the basis for detailed characterization of their physical properties (disks, binarity, 
atmospheres, accretion, activity). The survey makes use of extremely deep optical and near- 
infrared imaging, combined with the 2MASS and Spitzer photometry catalogs, and follow-up 
spectroscopy using 8-m class telescopes, aiming to detect the photosphere of the objects, 
which is an essential prerequisite for an unbiased survey. Our observations are designed to 
reach limiting masses of 0.003 M Q , well below the Deuterium burning limit at 0.015 M Q , 
and in each region, we aim to cover at least ~ lOOOarcmin 2 . 



In our first paper (IScholz et al.ll2009l ). we presented the survey of the NGC 1333 cluster 
in Perseus, which was surveyed complete down to 0.004-0.008 M Q for Ay up to 5-10 mag. We 
identified 12 new substellar members in this cluster based on photometry and spectroscopy, 
and found a tentative minimum mass cut off to the IMF at 0.012-0.020 M . 

In this second paper, we report on the surv ey of the p Oph star forming; cloud complex. 
p Oph is one of the closest (d = 125±25 pc, Ide Geus et al.l Il989l ) regions of active star 
formation. In contrast with the first SONYC field, NGC 1333, the p Oph cluster is not 
as compact and exhibits extremely high and variable levels of extinction. The main cloud, 
L168 8, is a dense molecular core, with visual extinction up to 50-100 mag (IWilking fc Lada 
19831 ). hosting an embedded infrared cluster of around 200 stars, inferred to have a median 
age o f 0.3 Myr, and surro unded by multiple clusters of young stars with a median age of 2.1 
Myr (IWilking et al.ll2005[ and references therein). 



1 Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observa- 
tory of Japan. 
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To date, a population of 17 BD candidates in p Oph, with spectral M6 o r later , have 
been spectroscopical ly confirmed by IWilking et al.l (119991 ). iLuhman &: Riekd (119991 ). and 
Cushing et all (hoooh An ISOCAM survey of very young stars in p Oph was presented in 
Bontemps et al.l (120011) , from which a sele ction of candidate BDs with disks was discussed in 
Natta et al.l (booj l. For a number of BDs, iJayawardhana et al I (I2002|) searched for ac cretion 
signatures with high resolution optical spectroscopy, and IJayawardhana et al.l (120031 ) found 
a disk fraction of ~67% for p Oph from L'-band (3.8 pm) e xcess. A survey of accretion rate 
measurements for the full ISOCAM survey was presented in lNatta et al.l ( 120061). For a recent 
summ ary of the star formation and brown dwarf population in p Oph, see IWilking et al. 
(120081 ). A few brown dwarfs were found and spectroscopically confirmed at about 1 degre e 
away from the center of the L1688 cloud (IJayawardhana fc Ivanovll2006l ; lAllers et al.l 120071 ) . 



This paper is structured as follows. Observations and data reduction of imaging and 
spectroscopy are described in Sect. [2j The photometric selection of very low mass objects is 
described in Sect. [3j and spectroscopy is presented in Sect. HI The results are discussed in 
Sect. [5j and conclusions are presented in Sect. [6j 



2. Multi-band photometry 
2.1. Optical imaging 



We obta ined optical images w ith the Subaru Prime Focus Camera (Suprime-Cam) wide 
field imager (IMiyazaki et al.l 120021 ) in the SDSS i' filter on August 2, 2008, with a typical 
range in seeing of 0.50-0.60". Suprime-Cam is a mosaic camera utilizing 10 CCDs arranged 
in a 5 x 2 pattern, giving a total field of view of 34' x 27', with a spatial resolution of 0.20" 
per pixel. The p Oph cluster was observed in a single field with 15 individual exposures of 
120 sec in a 10-point dither pattern. In total this gives an integration time of 1800 sec. 

We performed image reduction separate ly for each individual chip, using th e Suprime- 
Cam reduction software package SDFRED (jYagi et al.l 120021 ; lOuchi et al.l 120041 ) . This in- 



cludes overscan subtraction, flat-fielding, distortion correction, sky subtraction, bad pixel 
and AG probe m asking, and finally image combination. The 2MASS point source catalog 
(jCutri et al.ll2003l ) was used as a reference for the calibration of the world coordinate system, 
using the msctpeak routine from the IRAF package MSCRElfy. The typical fitting residuals 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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were of order 0.15". 

The objects were id entified using the Source Extractor (SExtractor) software package 



(IBertin Sz Arnoutslll996l ). We require an object to have at least 5 pixels with flux above the 
3cr detection limit to be extracted. The automatic aperture fitting routine in SExtractor 
was used to calculate the flux of each source. We chose conservative rejection criteria, 
to minimize the number of spurious detections in the photometry database, including the 
rejection of objects within 25 pixels of the edge of the image, elongated objects (a/b > 1.2), 
and saturated objects. The results from the detection and rejection algorithm were verified 
by visual inspection, to minimize the number of spurious detections and missed point- sources. 
The final optical catalog has 5791 objects. 

The chip-to-chip zeropoint offsets were derived from the median fluxes of domeflat im- 
ages. The absolute zeropoint for the mosaic was derived from observations of the SA112 



standard field, which contains SDSS secondary standards (ISmith et al.ll2002l ). These stan- 
dard observations were taken at the same airmass and with the telescope defocused to avoid 
saturation. Four standard stars were available in this field and photometry was obtained 
using a fixed 30-pixel radius aperture extraction. For the i'-band, we derive a mean zeropoint 
of 32.588 ± 0.076 mag, and a completeness limit of 24.15 ± 0.30 mag, see Fig. [TJ 



2.2. Near-infrared imaging 

The Multi-Object Infrared Camera and Spectrograph (MOIRCS) (jSuzuki et al. 2008 ) 



mounted on the Subaru telescope was used to observe p Oph in the J- (1.26 pm) and K s - 
band (2.14 pm). The data was obtained June 22 to June 24, 2007. MOIRCS uses two 
detector arrays which provide a field of view of 4' by 7' in each pointing. We observed 24 
MOIRCS fields with an overlap of 30" between adjace nt fields, encompa ssing; the embedded 



subclusters p Oph A, B, E, and F in the L1688 cloud (lLoren et al.lll990l ). The total imaging 
area coverage was 31.5' by 26'. Each field was covered in a six-point dither pattern with 
individual exposure times of 100 s in J-, and 50 s in K s -band. This yields a total integration 
time of 600 s in J- and 300 s in K s -band. Significant contamination by a bright nearby source 
was present in 15 individual images. 



The reduction was carried out with the SIMPLE-MOIRCS packagqj. First, we corrected 
for the offset in sensitivity between the two detector arrays. The individual frames are then 
flatfield corrected for flatfields and co-added by median. For the source detection we used the 



3 http://subarutelescope.org/Observing/Instruments/MOIRCS/imaginformation.htnil 




Fig. 1. — i'-band histogram of the objects in our photometric catalog with i'-band data. The 
peak at 24.15 mag indicates the completeness limit of the survey. The faintest objects in the 
survey are found at i > 26 mag. 
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Fig. 2. — J-band and K s -band histogram of the objects in our photometric catalog with J- 
and K s -band data. The peaks at 20.6 and 17.8 mag indicates the completeness limit of the 
survey in J- and K s -band respectively. 
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SExtractor, requiring at least 5 pixels with flux above the 3.5 a detection limit. Saturated 
and extended sources were rejected from further c onsideration. We calib rated the coordinate 



system against the 2MASS point source catalog (ISkrutskie et al.ll2006l ) for J-band and K s - 
band separately, with typical fitting residuals of 0.11". Duplicate sources from adjacent fields 
were removed. We cross-correlated the sources from the J-band catalog with the K s -band 
catalog, requiring a separation of < 2" to be a match. Sources that are only detected in one 
band, either J or K s , were rejected. The final near-infrared catalog contains 1571 objects. 

The fluxes for these sources were extracted using a fixed 6 pixel radius aperture for all 
fields. To correct for the variable seeing, we derived an aperture correction factor for each 
image and applied it to the fluxes. Based on the 2MASS point source catalog, we derived the 
absolute zero point for each individual field. For 19 fields no 2MASS sources were covered; 
for these fields a zero point is derived based on the average of zero points in other fields at 
similar airmass. The completeness limit of the observations was found to be 20.6±0.30 in 
J-band and 17.8±0.30 in K s -band, see Fig. [2j 



2.3. Cross correlation of optical and NIR catalogs 

The sources from the i'-band catalog were cross-correlated by coordinates with the 
sources from the JK S catalog, requiring a separation of < 2" to be a match. The resulting 
final catalog of sources with detections in all three bands (hereafter "iJK s catalog" ) contains 
504 sources. 

The spatial distribution of these 504 sources is shown in Fig. [3] in comparison with the 
visual extinction map given by the COMPLETE project^. As seen in the figure, we have 
chosen the pointings of our observations to cover the high extinction contours of the LI 688 
core; however, parts of the region surveyed by Suprime-Cam were not covered in MOIRCS. 

The majority of the iJK s catalog objects are located outside the region with the highest 
Ay, at a typical Ay range of 5-15. At higher Ay, toward the center of p Oph, the lack of 
sources corresponds to a lack of i'-band detections due to the high Ay. At lower Ay, many 
objects are saturated in JK S . As shown in Fig. [3j the majority of the BD candidates are 
found in several distinct clumps at regions with low to moderate cloud extinction, Ay = 
5-15, within the region covered by both the Suprime-Cam and MOIRCS imaging. Three of 
these regions are selected for multi-object follow-up spectroscopy, see Sect. HJ which showed 
that this apparent clustering is dominated by background objects. Thus, the distribution 



4 http: / /www.cfa.harvard.edu/COMPLETE 
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Fig. 3. — Spatial distribution of sources in p Oph. Contours are A v = 5, 10, 15, 20, 25, 30, 
as derived from 2MASS by the COMPLETE project. Dotted line: i-band imaging coverage; 
dashed line: J and K s -band imaging coverage; dash-dot line: MOIRCS spectroscopy fields. 
The locations of p Oph subclusters, as defined by the peaks of the DCO + dense cores from 



Loren et al.l (jl990l ). are indicated with letters. Top panel: i-band indicated +, J and K s - 
band indicated with x, iJK s catalog sources indicated with diamonds. Bottom panel: 
iJK s catalog sources indicated with +; iJK s catalog selected BD candidates with follow-up 
MOIRCS spectroscopy indicated with x (Sect.HJ); low mass YSOs spectroscopically confirmed 
in this paper (Sect. 14.21) indicated with red large squares. Candidate BD s selected from JK a + 
Spitzer photometry indicated with circles. Previo usly known BDs from lWilking et al.l (120081 ) 
and class II / III sources from lNatta et all ( 2006 ) are indicated with small blue squares and 
diamonds respectively. 
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of the combined iJK s catalog sources traces substructure in the extinction, and not in the 
distribution of YSOs. 

We convert the dynamic range and completeness limit of the i'-band and J-band survey 
separately to mass s ensitivity limits, by comparing with model evolutionary tracks DUSTYOO 
(kllard et all boOlh and COND03 faaraffe et al.l l2003h . described further in Sect. O and 
Figs. H and El 

The majority of i'-band magnitudes are between 18 and 25.5 mag, with the completeness 
limit at 24.15 mag. p Oph has significant and strongly variable extinction, with a typical 
minimum Ay of 5 mag, ranging up to 15-3 mag near the ce nter of the cluster, and upwards 
of 30 mag in t he very cente r of the cluster fjRidge et al.ll2006f ). For an Ay of 5, corresponding 
to Ai = 3.45 ( lMathislll990l ). the i'-band range of 18-25.5 mag corresponds to masses ranging 
between 0.002-0.068 M , with the completeness limit in i'-band corresponding to a mass 
ranging between 0.004 to 0.005 M based on the 1 Myr COND03 and DUSTYOO models. 
At a moderate Ay = 15, the i'-band survey range corresponds to masses from 0.047 to more 
than 0.1 M (outside model range), with the completeness limit corresponding to a mass 
sensitivity limit of 0.09 to 0.1 M . 

The J-band survey was sensitive to magnit udes between 15 and 21.3, with a complete- 
ness limit of 20.6 mag. For Ay = 5 (Aj = 3, lMathislll990l ) the 15-21.3 magnitude range 
corresponds to masses in the range of 0.001-0.034 M , with the completeness limit corre- 
sponds to masses of ~ 0.001 - 0.003 M . For an Ay of 15, the survey range corresponds to 
masses of 0.005-0.1 M , with the completeness limit corresponding to masses of 0.006-0.007 
M . 

In summary, the sensitivity of the i'-band and J-band survey were such that they were 
complete down to below the substellar regime, for low to moderate Ay, while sensitivity in 
the upper mass range was different between i'-band and J-band, due to J-band saturating 
at lower masses. However, particularly at large Ay, the survey was still not deep enough in 
i'-band to probe down into the substellar regime, and below. 



2.4. Cross correlation of optical and 2MASS NIR catalogs 



To increase the dynamic range of the survey in the n ear-infrared, we retrie ved additional 
JK S photometry from the 2MASS point-source catalog (ISkrutskie et al.l 120061 ) . We selected 
2MASS sources located within the spatial coverage of our i'-band survey. The objects were 
further required to not be extended with semi-major axis > 10", nor fall within the elliptical 
profile of such an extended source (xflag = 0), and not to be associated with known solar 
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system objects (a-flag = 0) and detection quality flags in J and K s -band of 'A'. The resulting 
2MASS catalog was cross-correlated by coordinates with the i'-band catalog sources. We 
require a separation of < 2" to be a match. The final i' + 2MASS J-, and K s -band catalog 
contains 264 sources. The majority of these sources were saturated in the MOIRCS data, as 
illustrated by Fig. HI 



2.5. Spitzer photometry 

We retrieved the Spitzer IRAC and MIPS photometry of p Oph from the Legacy Pro- 
gram data archive available at the Spitzer Science Center, using the "High reliability catalog" 
(HREL) created by the "Cores to Disks" (c2d) Legacy team, made available at the SSC web- 
sit^E The photometry from this catalog is used in Sect. [372] to search for candidate substellar 
objects with disks. 

The Spitzer catalog was cross-correlated with the MOIRCS J+K s -band catalog as well 
as the 2MASS catalog (see selection in Sect. 12.41) . Here we require a separation of < 5" for 
a match, and the final catalog is restricted to the spatial area of the Suprime-Cam i'-band 
coverage (246.415 < RA < 247.085, -24.77 < DEC < -24.275). The resulting catalog contains 
952 matches between MOIRCS and Spitzer, and 1166 matches between Spitzer and 2MASS. 



3. Selection of candidate low mass sources in p Oph 



3.1. iJK s candidate selection 



We initially selected low mass and BD candidates based solely on the (i', i'-J) colour- 
magnitude diagram. Fig. [4] shows the colour-magnitude diagram constructed from the 
Suprime-Cam i'-band catalog and bo th the MOIRCS a nd the 2MASS J-, + K s -band cata - 
log. For comparison, the DUSTY00 lAllard et alihoOlh and COND03 faaraffe et alil2003h 
atmosphere model evolutionary tracks for an age of 1 Myr are overplotted, adjusted to the 
distance of p Oph of 125 pc. We converted the Cousins I-band to the Sloan i'-band using the 
transformation given by iJordi et ajj (120061 ). The BD candidate selection cut-off in (i'-J) is 
shown in Fig. [4] and was determined as a linear approximation of the 1 Myr COND03 track. 
Any object redder than this line is considered a candidate. This selection gives us 309 BD 
candidates based on our iJK s catalog and an additional 228 sources based on the Suprime- 



5 http://data. spitzer. caltech.edu/popular/c2d/20071101_enhanced_vl/oph/catalogs/ 
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Cam i'-band + 2MASS JK S catalog. We obtained follow-up spectroscopy for a selection of 
the iJK s catalog based BD candidates, as presented in Sect. HI 

The intrinsic J-K s color of late-type objects is mostly independent of spectral type 
and luminosity and can be used to estimate the extinction for suspected cluster members. 
We adopt here an intrinsic color of J-K s = 1.0 mag, which is consistent, within ± 0.2 
ma g, with the empirical values for main sequence dwarfs and giants with spectral type 



M f Bessell fc Brett 



19881 ). with the predictions from the evolutionary tracks COND03 (for 



0.003-0.1 M ), BCAH98 (for 0.03-0.5 M ), and DUSTY00 (for 0.006-0.1 M ? ), and with the 



color s of BDs in the 3 Myr old a Ori star forming region down to 0.013 M (ICaballero et al. 



20071 ). For young objects with masses be low the Deuterium burning limit, there are indica- 
tions that the J-K s color increases to 1.5 ( ILodieu et al.ll2006t ICaballero et al.ll2007l ); similarly 
the presence of a circumstellar disk can increase the J-K s color. 



To estimate the extinction we use the reddening law from iMathid (119901 ) with Ry = 4. 
Fig. [5] shows the (J, J— K s ) colour-magnitude diagram for the sample of 309+228 i'-band + 
MOIRCS and 2MASS catalog derived BD candidates. The reddening path is indicated in 
solid and dashed lines, and the intrinsic color indicated with a vertical dotted line. The 
lMyr COND03 evolutionary track is used for the mass estimate based on the de-reddened 
J-band magnitude. We obtain visual extinctions of Ay = 3-16 for the MOIRCS sources 
and A v = 2 - 34 for the 2MASS sources. 

The main errors in the extinction estimates are the uncertainties in the estimate of 
the intrinsic colors (accurate to within 0.2 mag in J-K s ) and the po ssible excess emission 



due to disks and accretion (< 1 mag in J-K s , see iMeyer et al.lll997f). par t icular ly for the 
sample selected based on Spitzer MIR excess in Sect. 13.21 IWilking et al.l (120081 ) lists 316 
L1688 association members known in January 2008, out of which 152 show infrared excess 
indicative of the presence of a disk, (indicated with 'B2ex', 'B4ex', and 'Mex' in their Table 
1). This illustrates that a substantial fraction of the L1688 members may show enhanced 
K s -band flux due to excess emission from the dust in the disk. We conservatively estimate 
that our typical uncertainty in Ay is in the range of ±1 mag for most sources, but may be 
up to 5 mag if disk excess is present. We note that brown dwarfs with disks typically have 
only little colour excess in the near-infrared. 



The estimated range of visual extinction is in good agreement with the p Oph extinction 
map as derived from 2MASS imaging observations by the COMPLETE team. 
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Fig. 4. — Colour- magnitude diagram in (i', i'-J). Sources from the Suprime-Cam and 
MOIRCS iJK s catalog are indicated with light "+" symbols, while BD candidates based 
on iJK s are indicated in dark "+" symbols. Sources from the Suprime-Cam and 2MASS cat- 
alogs are indicated with purple "+" symbols. Spectroscopy follow-up targets are indicated 
with x, two confirmed ver y low mass YSOs with large red squares. Previously known BDs 
from IWilking et al.l (120081 ) are indicated with small blue squares, class II / III sources from 



NattaetaL 



( 120061 ) are indicated with small blue diamonds. A tmosphere mo del isochrones 



for ag e 1 Myr are overplotted for DUSTYOO flAUard et al.ll200ll ) and COND03 (IBaraffe et al. 



20031 ). The solid line denotes the BD candidate selection cut-off. 
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Fig. 5. — Colour-magnitude diagram in (J, J-K s ), constructed from the MOIRCS and 
2MASS photometry, for the BD candidates selected from the iJK s catalog and the i-band + 
2MASS catalog. Solid and dashed lines indicate visual extinction of 10 and 20 magnitude 
from the assumed intrinsic J-K s of 1 (vertical dotted line). "+" and diamond symbols indi- 
cate the MOIRCS and 2MASS photometric points, circles are dereddened in J. The two very 
low mass candidates confirmed by MOIRCS spectroscopy are included with red symbols. 
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3.2. JK S + Spitzer candidate selection 

As the i'-band observations were not complete in the central parts of L1688 with high 
Ay, we selected a second sample of candidate substellar mass objects with disks based on 
the mid-infrared photometry from the Spitzer Telescope, restricted to the spatial area of the 
i'-band coverage (see Sect. 12. 5p . Fig. |H] shows a colour-colour diagram of the 4 IRAC bands 
for the MOIRCS and 2 MASS sources. For 11 MOIRCS sources, and 72 2MASS sources, the 
IRAC colours are consistent or redder than sources with a circumstellar disk. The presence 
of a circumstellar disk is taken as a confirmation of youth and membership of the p Oph 
cluster. 

Fig. [7] shows the (J, J-K s ) colour-magnitude diagram for the sample of 11 MOIRCS and 
72 2MASS sources with excesses in the Spitzer data. Based on the CMD and the DUSTYOO 
atmospheric model, we find 10 MOIRCS BD candidates and 17 2MASS BD candidates with 
M < 0.08 M , assuming an intrinsic brown dwarf photosphere J-K s color of 1 magnitude. 
These 10 + 17 candidates are listed in Table [TJ Among this selection of 27 BD candidates, 16 
are new, while 11 are previously spectroscopically confirmed BDs. If we assume an additional 
1 magnitude J-K s excess due to the presence of a circumstellar disk (intrinsic J-K s = 2), 
then the number of candidates increases to 11 MOIRCS and 39 2MASS sources. 
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Fig. 6. — Colour-colour diagram in (3.6-4.5, 5.8-8.0), constructed from Spitzer IRAC pho- 
tometry, for 2MASS (crosses) and MOIRCS (diamonds) sources with a Spitzer match. Th e 
dashed line denotes the area where class II objects are located, based on lAllen et al.l ( 12004] ) . 
All sources selected as having significant NIR excess are indicated with in black. 



Table 1. Likely substellar members with disks within i'-band surveyed area within p Oph from Spitzer + JK S 

photometry. 



# 




RA 




Dec 


J 


K s 


IRAC1 


IRAC2 


IRAC3 


IRAC4 




JK S phot. 


Identifiers 


References' 1 




(J2000) 


(J2000) 


(mag) 


(mag) 


(mjy) 


(mjy) 


(mjy) 


(mJy) 


(mag) 






BD candidacy 


1 


16 


26 18.58 


-24 


29 51.4 


17.36 


13.61 


3.93 


4.08 


4.58 


5.44 


14.5 


MOIRCS 


BKLT J162618-242951 


- 


2 


16 


26 22.28 


-24 


24 07.1 


16.70 


13.94 


3.10 


4.26 


5.53 


9.21 


9.3 


MOIRCS 


GY 11 


WGM99,CTK00,N02 


3 


16 


26 36.84 


-24 


19 00.1 


17.01 


14.55 


1.10 


2.05 


3.70 


3.50 


7.7 


MOIRCS 


BKLT J162636-241902 


- 


4 


16 


26 56.37 


-24 


41 20.5 


18.70 


14.77 


1.01 


0.99 


1.06 


1.07 


15.5 


MOIRCS 




- 


5 


16 


26 58.66 


-24 


24 55.6 


19.64 


15.30 


3.72 


5.17 


5.65 


5.84 


17.6 


MOIRCS 


AOC J162658.65-242455.5 


- 


6 


16 


26 59.04 


-24 


35 56.9 


16.51 


12.21 


14.20 


13.70 


11.60 


10.40 


17.5 


MOIRCS 


GY 172 


- 


7 


16 


27 32.15 


-24 


29 43.6 


18.47 


13.41 


8.28 


9.14 


9.16 


8.53 


21.6 


MOIRCS 


GY 287 




8 


16 


27 38.95 


-24 


40 20.7 


16.54 


13.36 


15.80 


26.70 


38.20 


53.20 


11.5 


MOIRCS 


GY 312 


- 


9 


16 


27 41.80 


-24 


42 34.7 


21.32 


14.97 


1.64 


2.03 


2.29 


2.43 


28.3 


MOIRCS 


- 


- 


10 


16 


27 48.24 


-24 


42 25.8 


20.76 


15.21 


1.21 


2.21 


3.64 


14.10 


24.1 


MOIRCS 


AOC J162748. 24-244225. 6 


- 


11 


16 


25 42.54 


-24 


43 51.1 


15.61 


14.03 


0.99 


0.66 


0.47 


0.39 


3.0 


2MASS 


BKLT J162542-244350 


- 


12 


16 


26 04.58 


-24 


17 51.5 


15.79 


12.19 


10.20 


9.89 


8.98 


11.30 


13.8 


2MASS 


BKLT J162604-241753 


- 


13 


16 


26 18.82 


-24 


26 10.5 


14.84 


12.14 


10.50 


11.50 


13.90 


20.50 


9.0 


2MASS 


CRBR 2317.3-1925 


WGM99 i- 


11 


16 


26 21.53 


-24 


26 01.0 


12.57 


10.92 


25.00 


20.30 


19.30 


19.10 


3. 1 


2MASS 


GY 5 


WGM99,N02 c 


15 


16 


26 21.90 


-24 


44 39.8 


12.34 


10.86 


28.80 


26.60 


23.80 


26.40 


2.6 


2MASS 


GY 3 


N02,W05 


16 


16 


26 23.81 


-24 


18 29.0 


16.07 


13.58 


2.72 


2.27 


1.79 


1.62 


7.9 


2MASS 


CRBR 2322.3-1143 


CTK00 


17 


16 


26 27.81 


-24 


26 41.8 


14.26 


12.09 


9.52 


8.74 


7.78 


7.38 


6.2 


2MASS 


GY 37 


WGM99,W05 


18 


16 


26 37.81 


-24 


39 03.2 


14.98 


12.85 


4.47 


3.92 


3.35 


3.21 


6.0 


2MASS 


GY 80 




19 


16 


26 51.28 


-24 


32 42.0 


15.30 


13.89 


1.55 


1.31 


1.34 


1.24 


2.2 


2MASS 


GY 141 


LLR97.CTK00 


20 


16 


27 05.98 


-24 


28 36.3 


16.64 


12.97 


5.43 


5.16 


5.32 


5.37 


14.1 


2MASS 


GY 202 


WGM99,LR99,CTK00 


21 


16 


27 26.22 


-24 


19 23.0 


16.40 


12.93 


6.20 


6.29 


6.48 


6.83 


13.1 


2MASS 


BKLT J162726-241925 




22 


16 


27 26.58 


-24 


25 54.4 


13.00 


11.84 


10.30 


10.20 


9.78 


12.00 


0.9 


2MASS 


GY 264 


W05 


23 


16 


27 37.42 


-24 


17 54.9 


14.15 


11.95 


11.20 


10.20 


10.20 


11.70 


6.4 


2MASS 


ISO-Oph 160 


N02 


24 


16 


27 38.95 


-24 


40 20.7 


16.54 


12.29 


15.80 


26.70 


38.20 


53.20 


17.2 


2MASS 


ISO-Oph 165 




25 


16 


27 40.84 


-24 


29 00.7 


14.66 


13.10 


2.78 


2.28 


2.22 


2.64 


2.9 


2MASS 


GY 320 




26 


16 


27 46.29 


-24 


31 41.2 


13.83 


11.32 


17.60 


15.70 


14.80 


18.60 


8.0 


2MASS 


GY 350 


N02 


27 


16 


28 17.74 


-24 


17 05.0 


14.57 


12.31 


7.40 


6.17 


5.20 


4.29 


6.6 


2MASS 


BKLT J162817-241706 





a References for pa pers spectroscopically confirming BP candidacy: WGM 99: IWilking et"al"l il999l). CTK 00: ICushing et"al] feood), N0 2: iNatta et all J2002ll. W05: 
IWilking etall 11200511. References for identifier s: GY: iGreene fc Young! jl992f) : CRBR: IComeron et al I h993h ; BKLT: iBarsonv et al.l Jl997h ; ISO-Oph: iBontemps etall 
fcOQll) ; AOG lAlves i de Oliveira fc Casalil feOQSl l 
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For this sample of JK S + Spitzer excess sources, estimates of the extinction are derived 
as described in Sect. 13.11 The visual extinction obtained for this sample is Ay = — 30, and 
is given in Table [TJ 



4. Multi-object spectroscopic follow-up of candidates 

4.1. Observations and data reduction 

We used MOIRCS to carry out multi-object spectroscopy for 59 BD candidates selected 
from the i', J, K s photometry (see Sect. [3]), comprised of 58 new candidates and 1 previously 
known BD candidate, GY 84. The targets were covered in 3 fields, centered on the areas 
with highest density of candidates, indicated in Fig. [3J The three fields are numbered by 
increasing RA as field 1, 2, and 3. Two of these candidates are located only 1.5" apart 
from each other in a, so that they can be cover ed in one sli t. In a ddition we observed one 



confirmed young M dwarf, GY84, published by iNatta et all (120061 ). which falls into one of 
our fields. 

Pre-imaging for the MOS masks in K s -band was obtained in April 2009. We used slits 
which are 0.9" wide and 9-12" long. The spectroscopy run was carried out in two nights on 
May 30-31, 2009 using the grisms HK500, which covers the H- and K s -band. The second 
half of the first night was lost due to technical problems, as a result the exposure time for 
field 1 is lower than intended. The total on-source time was 30, 140, and 150 min for fields 
1-3, split in shorter exposures of 5 or 10 min. The typical seeing during these nights was 
0.5", as determined in the focusing procedure. The average airmass for the three fields is 
1.8, 1.5, and 1.5. 

Before and after the p Oph masks, we observed AO stars through one of the science 
masks for flux calibration and extinction correction; these exposures required de-focusing to 
prevent saturation. In total two AO stars were observed in night 1 and seven in night 2, 
at various airmasses. For each mask, we obtained series of domeflats with lamp on and off 
for calibration purposes. Between the single exposures, the telescope was moved so that the 
targets shifted by 2.5" along the slit (nodding), to facilitate sky subtraction. 

The MOS data was reduced following mostly the recipes outlined in Scholz et al. (2009). 
After subtracting the nodded exposures to remove the background, the images were divided 
by a normalized flatfield. Frames from the same nodding position were co-added using the 
median and a sigclip routine to remove cosmic ray events. This gives us two final frames 
for each chip and mask, from which the spectra were extracted using apall in IRAF. All 
candidates are well-detected. The wavelength solution was derived from the OH lines in the 
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Fig. 7. — Colour- magnitude diagram in (J, J-K s ), constructed from the MOIRCS and 
2MASS photometry, for the sources with Spitzer NIR excess. Solid and dashed lines in- 
dicate visual extinction of 10 and 20 magnitude from the assumed intrinsic J-K s of 1. "+" 
and diamond symbols indicate the MOIRCS and 2MASS photometric points, circles are 
dereddened in J. 
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unreduced frames, for each object separately. The typical RMS of the fit was 2-3 A, well 
below the resolution. All spectra binned to 40 A per pixel, corresponding to a resolution of 
500 at 2.0 /im. The two nodded spectra for the same object were co-added. The standard 
star spectra were treated in exactly the same way as the science frames. 

After background subtraction the science frames still show 5-10% residuals from the OH 
lines. The reason is that p Oph was observed at relatively high airmasses, so that consecutive 
exposures differ in airmass by as much as 0.1. However, the extraction with apall includes a 
background fit perpendicular to the dispersion direction which reliably removes the residuals. 
We made sure that the final science spectra are not significantly affected by OH lines. 

The corrections for the broad telluric absorption features between the near-infrared 
bands and for the instrument response was done in one step, using the AO s tandard stars . 



The AO standard star spectra were divided by a library spectrum of a AO star (lPickleslll998l ). 
obtained from the ESO website^]. The resulting calibration spectra show the combined effects 
of telluric absorption and response. Two out of the nine spectra differ significantly from the 
average; in one case probably due to effects of saturation, the other one has an usually large 
J-K s colour of 0.2 (obtained from 2MASS) which might indicate the presence of a dust shell. 
Both were not considered for the calibration. 

Most of the differences in the remaining 7 AO spectra are caused by the variable airmass. 
For each wavelength bin we fit the fluxes linearly as a function of airmass and extrapolated to 
airmass of 1.0. While the extinction coefficients are consistent with the standard extinction 
law in the wavelength domain 1.45-1.8 pm and 2.05-2.3 /im, they vary strongly and irregularly 
in the water absorption features outside these bands, which cannot be reliably calibrated. 
The spectrum for airmass 1.0 was smoothed and fit with a second order polynomial over 
the H- and K s -band. The result is our final calibration spectrum. All science spectra were 
extinction corrected using the coefficients derived from the standard spectra and the average 
airmass, and then divided by the calibration spectrum. 



4.2. Spectral analysis 



Our goal is to identify y oung sources with spectral types later than M5. As outlined in 
detail in lScholz et al.l (120091 ). these objects show a characteristic spect ral shape in the near - 
infrared. In particular, their spectra have a clear peak in the H-band (jCushing et al.ll2005l ). 
This feature is caused by water absorption on both sides of the H-band. The depth of these 



http: / /www. eso.org/ sci/facilities /paranal/instruments /isaac / tools /lib / 
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features depends strongly on effective temperature. While the H-band peak appears round 
in old field dwarfs, it is unambiguous triangu lar in young, low-gravity sources, likely du e to 
the effects of Collisional Induced Absorption ( IKirkpatrick et al.ll2006l ; IWitte et al.l 120091 ). In 
addition, young brown dwarfs are expected to have flat or increasing K s -band spectra with 
CO absorption bands at A > 2.3 pm. 

In a first step we visually inspect our spectra to search for the characteristic signature 
of young brown d warfs . Three ob j ec t s exhibit the features described above, including the 
object GY 84 from lNatta et al.l ( 120061 ). these objects are listed in Tableland their spectra are 
shown in Fig. [HJ For six objects the signal-to-noise ratio is too low to conclusively determine 
their nature. The remaining 50 objects have smooth spectra with mostly decreasing K s -band 
slope. 

In a second step we estimate the e ffective temperatur es of the BD candidates by fitting 
m odel spectra f r om th e DUSTY series (lAllard et al.ll200ll ). following the procedure outlined 



m iScholz et al.l (120091 ) . We used the extinction as deter mined from the J- K s colour (see 
Sect. HO}. For GY84 the value of A v = 14 published bv iNatta etaJ tood ) is used. The 
results are listed in Table [21 We find best-fitting effective temperatures of 2500 and 3100 K 
for our two candidates, here named SONYC-RhoOph-1, and SONYC-RhoOph-2, and a T fj 
= 3400 K for GY84, corresponding to spectral types M9, M5, and M3 respectively. SONYC- 
RhoOph-1 is confirmed as a newly found brown dwarf. Based on COND03 and DUSTY00 
evolutionary tracks for lMyr, it has a mass estimate of 0.015-0.02 M , marking it as one 
of the lowest mass objects in p Oph. The temperatures we derive for G Y 84 and SONYC- 
RhoOph-2 are too high to classify them as brown dwarfs. For comparison INatta et al.l (120061 ) 
list T e fj = 2900 for GY84, derived by converting the J-band magnitude to luminosity and 
comparing with theoretical evolutionary tracks. The model spectrum with 2900 K however 
clearly does not match the observed spectrum. 

The remaining sources, which are not brown dwarfs and have sufficient signal-to-noise 
ratio, exhibit smooth, featureless spectra, which makes it difficult to constrain their nature. 
The overwhelming majority of these objects are expected to be background late-type objects 



Table 2. Probable low mass and substellar mass members of p Oph, with MOIRCS 

spectroscopy follow-up. 



# 


RA (J2000) 


Dec (J2000) 


i' (mag) 


J (mag) 


K s (mag) 


T cff (K) 




Notes 


1 


16 26 56.33 


-24 42 37.8 


21.24 


17.68 


15.53 


2500 


5 


SONYC-RhoOph-1 


2 


16 26 57.36 


-24 42 18.8 


22.53 


18.94 


15.73 


3100 


10 


SONYC-RhoOph-2 


3 


16 26 38.82 


-24 23 24.7 


21.10 


15.10 


sat 


3400 


14 


GY 84 
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Fig. 8. — MOIRCS spectra (solid) of candidate substellar objects in p Oph, with best fit 
reddened model (dashed) with an T eff of 2500 K, 3100 K, and 3400 K; offsets were applied 
for clarity. 
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giant or dwarf stars. Some of them could be embedded young stellar objects with spectral 
types earlier than M, although we do not have any evidence for youth. We note that this 
sampl e contains 14 objects^from the BKLT catalog, a near-infrared survey by Barsony et al. 



( 119971 ) without spectroscopic follow-up. Four examples of these featureless spectra are shown 
in Fig. El 



4.3. Comparison with other studies 



The iJK s selected cata log of candidate brow n dwarfs is compared with the list of 24 
candidate BDs published in IWilking et al.l (120081 hereafter W08). 21 objects fall within the 
spatial coverage of both the i'-band and the J+K s -band imaging; 3 objects are not covered 
in i'- and/or J-band (GY141, Oph-160, Oph-193). Of the 21 covered, 12 are saturated 
(CRBR14, CRBR31, GY3, GY5, GY10, GY37, GY64, GY84, GY204, GY264, GY310, 
GY350, GY202), 3 are undetected in i-band (CRBR15, CRBR23, GY258), 1 is undetected 
in both i-band and J-band (GY31), and 1 is not properly detected in J-band due to strong 
distortion from a nearby bright source (see Sect. 12 .2 j) . 3 of the 21 are recovered in the iJK s 
catalog, GY11, GY201, and GY325, indicated in Fig. HI Of these 3, only 1 is recovered as a 
BD candidate (GY11) in the (i\ i'-J) CMD selection (Sect. EH), while the other 2 (GY201, 
GY325) are blue-ward of the selection cutoff. 



We also compare our candidate BDs to the iNatta et al.l (120061 ) sample of objects with 
evidence of a circumstellar disk, and measurements of accretion. Three sources are covered 
in the iJK s catalog, ISO033 (=GY11 in W08), IS0165, and IS0169b, all three are recovered 
as brown dwarf candidates. 



We compare our study with the recent paper by lAlves de Oliveira et al.l (120101 ). hereafter 
referred to as A10. Both studies use JHK S near-infrared photometry; in this paper, i-band 
photometry is added for the entire survey area, while A10 added archive i-band and z-band 
photometry for a subset of their JHK S sources. Both studies select candidate BDs from 
colour-magnitude and colour-colour diagrams based on predictions by evolutionary models. 
Two of the six new BDs found by A10 are covered in our iJK s survey area, but were not 
recovered in our iJK s catalog, due to saturation in J-band in one case and extended emission 
interference in i'-band in the other case. Our one new spectroscopically confirmed BD (RA 
246.734 -24.711) appears to have been covered by A10 (based on their Fig. 4), but is not 



7 BKLT J162640-242046, J162646-242155, J162654-244254, J162655-244242, J162656-244338, J162658- 
244342, J162703-244400, J162704-244422, J162708-244227, J162743-242822, J162745-242952, J162748- 
242542, J162750-242544, J162759-242912 
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Fig. 9. — Example MOIRCS spectra (solid) of sources in p Oph which are not substellar 
candidates, with overplotted reddened blackbodies (dashed) with T e fr ranging from 3000- 
4000 K; offsets were applied for clarity. The objects #1 and #4 in this plot are known in 
the literature as BKLT J162759-242912 and BKLT J162708-244227 faarsonv et alill997h . 
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reported as a candidate Rho Oph member (A10, Table 4). 



5. Discussion: brown dwarf population in p Oph 



In the following section we use our survey results in combination with literature data 
to put limits on the total number of substellar objects in the p Oph star forming region. We 
would like to preface this analysis with a cautionary note: p Oph with its high and strongly 
varying extinction is a difficult field to search for brown dwarfs. All existing surveys have 
to be considered incomplete. The three key requirements to achieve a full census of BDs in 
this area are a) observations in multiple bands, b) combination of different survey strategies 
(e.g., broadband vs. narrowband surveys), and c) comprehensive spectroscopic follow-up. 

Our SONYC survey provides new constraints on the substellar population in p Oph. We 
have used two different strategies. In a first step, we identify objects with iJK s photometry. 
As it turns out, this selection is heavily affected by background contamination and is not 
effective in the high extinction areas. It mostly covers the edges of the dense cloud core. In 
a second step we look for Spitzer/IRAC colour excess due to a disk in objects from our JK S 
database. This sample should be dominated by young objects in p Oph, but it is not clear 
what fraction is substellar. Since we do not require an optical counterpart, this selection also 
traces the dense cloud areas. Thus we provide two complementary approaches to identify 
brown dwarfs. 

Our iJK s survey is certainly not complete for substellar objects. In the high extinction 
areas the i-band data is not deep enough (although deeper than all other previously published 
optical surveys in this region). In the low extinction areas, more massive brown dwarfs can 
be saturated in the J- and K s -band, as evidenced by the significant number of literature 
BDs missing in our catalog (Sect. 14. 3p . Our sample identified from Spitzer colour excess is 
missing objects without disks. 

The previous surveys in the literature are s ubject to various biases as well. For example, 
the 20 objects observed by lWilking et al.l (119991 ) have been selected from various near-infrared 
surveys which yield a muc h larger initial sampl e. Their colour criterion favours objects 'near 



Natta et al. 



the surface of the cloud' ( Wilking et al.lll999r) and will not find deeply embedded sources. 
The objects confirmed by 



(2002) as br own dwarfs in p Oph have been selected 
from the ISOCAM survey by lBontemps et al.l (120011 ) using a combined extinction (i.e. colour) 
and luminosity criterion. This approach will by definition not be able to find the objects 
without disks nor the ones with high extinctio n s. A simila r argument c an be made for the 
objects confirmed by IJayawardhana &: Ivanovl (120061 ) and lAUers et al.l (120071 ). which were 
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selected from Spitzer data. This simply illustrates that only the combination of multiple 
survey strategies can provide a more complete census. 

One of the key features of our survey is extensive follow-up spectroscopy for our primary 
candidates from the iJK s database. Out of 504 objects, 309 were selected with colours 
consistent with substellar or planetary masses; for 58 of these candidates, and 1 known 
YSO, we obtained spectra. Only three of these turn out to be bona fide young stellar objects, 
one is a newly discovered young brown dwarf. Assuming that we selected the spectroscopy 
candidates in an unbiased way, we expect that the total catalog contains not more than a 
handful (~ 5) of new substellar objects. 

The latter estimate is only a lower limit for the actual number of new BDs in our 
catalog. The spectroscopy fields were deliberately placed in regions with strong clustering. 
As the spectroscopy shows, most of the objects in these regions are background stars, i.e. 
the clustering reflects substructure in the extinction of the cloud. In areas with higher or 
lower extinction we expect to have less contamination by background objects because they 
are either not visible or do not pass our colour criterion. Taking this into account, we make 
a second estimate of the total number of BDs in the iJK s catalog based on the success rate 
and spatial coverage of the follow-up spectroscopy while we assume a spatially homogeneous 
population of brown dwarfs. Our three MOS fields cover 72 sq. arcmin out of ~ 616 sq. 
arcmin (0.171 sq. deg.) in the total i+JK s survey (~ 12%). We obtained spectra for 59 out 
of the 159 iJK s candidates (~ 37%) located within the area of the MOS fields and find one 
bona fide brown dwarf. Scaling to the total area of our survey yields an estimate of 24 brown 
dwarfs. We note this estimate is taken as an order of magnitude estimate, since it has a large 
uncertainty (1 a confidence interv al of 5-74, assuming bi nomial distribution), and moderate 



clustering in the YSO population ( jBontemps et al.ll200fl ) may further increase the estimate. 



We conclude that despite the low yield of our spectroscopy campaign our database might 
still contain a substantial number of new brown dwarfs. Spectroscopy covering the full field 
is required to clarify this. 

We demonstrate that the available Spitzer data contains more objects than previously 
known that might be BDs in p Oph. By combining our deep JK S data with the Spitzer C2D 
catalog, we select 10 candidates with mid-infrared colour excess and near-infrared colours 
indicative for a substellar mass, of which 1 is a previously spectroscopically confirmed BD. 
In addition, 17 objects were found with 2MASS JK S data in combination with Spitzer data, 
using the same approach, of which 10 objects were previously spectroscopically confirmed 
BDs. The 17 new candidates exhibit a wide range of visual extinctions from to 20 mag. In 
particular the high extinction objects constitute a parameter regime not covered sufficiently 
in previous surveys. Again, this sample requires full spectroscopy follow-up for a more 
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detailed characterisation. 



In Table [3] we provide a list of the spectroscopically confirmed brown dwarfs in the 
literature, including our new object, restricted to our i'-band survey area (246.415 < RA < 
247.085, -24.77 < DEC < -24.275). The main criteria for an object to appear in this list are 
a) spectral type later than M6 or alternatively effective temperature below 3000 K and b) 
evidence for youth (e.g., disk excess, spectral features indicative of low gravity). We include 
sources with multiple spectroscopically determined spectral type estimates, if at least one is 
M6 or later. To date, 21 objects fulfill these criteria. As outlined above in detail, this number 
is almost certainly a lower limit to the full number of BDs in p Oph. The table should be 
seen as a starting point in establishing a m ore complete census. We did not include the 
three brown dwarfs from lAUers et al.l (120071 ) as those are located outside our survey area. 



We no te the inclusion of the spectroscopically confirmed early T-dwarf found by lMarsh et al. 



(120101 ). for which evidence for youth is found in low gravity features. This may be the lowest 
mass object identified in p Op h so far, although its cluster membership is still i n question as 
Alves de Oliveira et al.l ( 20101 ) note that they find a different K s magnitude than lMarsh et al. 



(120101 ). based on which the distance estimate would place the object at 137 - 217 pc, behind 
the p Oph cloud. Finally, we note that 16 of the 21 objects in Table [3] are in agreement 
with lAlves de Oliveira et al.l ( I2010l . Table 7). We do not list the 5 of their BDs outside our 
survey area: ISO-Oph 193, CFHTWIR-Oph 4, 47, 57, and 106. We do list 5 BD candidates 
that they do not : the newly found SONYC-RhoOph-1 and MARSH10-4450, as well as the 
3 borderline cases GY 37, GY 59, and GY 84, that have multiple conflicting spectral type 
estimates in the literature. 



Table 3. Spectroscopically confirmed Brown Dwarfs within i'-band survey area of p Oph. 



Name(s) RA (J2000) Dec (J2000) Indicator(s) of youth SpT References' 1 



CRBR 2317.3-1925 / CRBR 14 


16 


26 


18 


82 


-24 


26 


10 


5 


CO bands+Av, mid-IR excess 


M7.5,M5.5,M7 


WGM99,LR99,N02 


GY 5 


16 


26 


21 


54 


-24 


26 


01 





CaH, TiO, CO bands+Av, mid-IR exc. 


M5.5,M7,M6 


W05,WGM99,N02 


GY 3 


16 


26 


22 


05 


-24 


44 


37 


5 


CaH, TiO, mid-IR excess 


M8,M7.5 


W05,N02 


GY 10 


16 


26 


22 


17 


-24 


23 


54 


4 


CO bands+Av 


M8.5.M6.5 


WGM99.LR99 


GY 11 


16 


26 


22 


28 


-24 


24 


09 


3 


H 2 abs., CO bands+Av 


M6.1,M6.5,M8.5 


CTK00,WGM99,N02 


CRBR 2322.3-1143 / CRBR 31 


16 


26 


23 


78 


-24 


18 


31 


4 


H 2 abs. 


M6.7 


CTK00 


GY 37 


16 


26 


27 


83 


-24 


26 


42 


6 


CaH, TiO 


M5,M6 


W05,WGM99 


GY 59 


16 


26 


31 


37 


-24 


25 


30 


3 


CO bands+Av 


M3.75,M6,M5 


W05,WGM99,LR99 


GY 64 


16 


26 


32 


56 


-24 


26 


36 


9 


H 2 abs., CO bands+Av 


M7.0.M8 


CTK00,WGM99 


GY 84 


16 


26 


38 


80 


-24 


23 


22 


7 


CO bands+Av, H2O abs. 


M6,M3 b 


WGM99, this work 


CFHTWIR-Oph 34 


16 


26 


39 


92 


-24 


22 


33 


6 


H 2 absorption 


M8.25 


A10 


GY 141 


16 


26 


51 


42 


-24 


32 


42 


7 


H a , H2O abs. 


M8.5,M8.0 


LLR97,CTK00 


SONYC-RhoOph-1 


16 


26 


56 


33 


-24 


42 


37 


8 


H2O absorption 


M9 b 


this work 


GY 202 


16 


27 


06 


00 


-24 


28 


37 


3 


H 2 abs., CO bands+Av 


M5.7,M7,M6.5 


CTK00,WGM99,LR99 


GY 204 


16 


27 


06 


58 


-24 


41 


17 


9 


CaH, TiO, mid-IR excess 


M5.5,M6 


W05,N02 


MARSH10-4450 C 


16 


27 


25 


35 


-24 


25 


37 


5 


CH4, H2O absorption 


early T 


M10 


GY 264 


16 


27 


26 


58 


-24 


25 


55 


1 


CaH, TiO 


M8 


W05 


ISO-Oph 160 


16 


27 


37 


30 


-24 


17 


56 


1 


mid-IR excess 


M6 


N02 


GY 310 


16 


27 


38 


67 


-24 


38 


38 


2 


CO bands+Av 


M8.5,M7,M6 


WGM99,LR99,N02 


CFHTWIR-Oph 96 


16 


27 


1Q 


84 


-24 


29 


00 


8 


H2O absorption 


M8.25 


A10 


GY 350 


16 


27 


16 


36 


-21 


31 


11 


6 


disk excess 


M6 


N02 



1-0 



a R 97: iLuhman et all Jl997h . W GM99: IWilking et aL il999h , LR99: iLuhman fc Riekel lll999h , CT K00: ICushing et al. ll2000f) , N02: iNatta et all [|2002l '), 
W05: IWilking et al.l \200d ), W0 8: IWilking et alj J2008h . A10: lAlves de Oliveira et al.l feoicl) , M10: iMarsh et all fcOlCll References for identifiers: GY: 
iGreene fc Yound l|l992f ). CRBR: IComeron et alj jl993h . ISO-Oph: iBontemps et al.1 |200lj i 

b Spectral type estimate based on T e ff - SpT scale from Table 8 of lLuhman et al . (2003) 



c Provisional name chosen in this paper; no entry for this source exists in SIMBAD 
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Our sample of candidates identified from Spitzer excess allows us to estimate an upper 
limit for the number of missing BDs in p Oph. Our MOIRCS JK S survey area covers 19 
of the 21 spectroscopically confirmed brown dwarfs from Table El In the same area are an 
additional 16 new BP cand idates with disks (see Tabled]). Assuming a disk fraction of ~ 67% 
( jJayawardhana et al.l 120031 ). the upper limit for BD candidates in this area is ~ 16 x 1.5 = 
24, and subsequently, an upper limit on total number of BDs would be 24 + 19 = 43, roughly 
twice the number already confirmed with spectroscopy within our JK s -survey area. In other 
words, the current census of confirmed BDs might still miss a few tens of substellar objects. 

It is premature to derive the substellar IMF in p Oph. However, based on the currently 
available data we can estimate the number of low-mass stars (0.1 < M < 1.0 M ) in relation 
to the number o f BDs (M < 0.1 Mm) , a ratio that has been derived for a number of star 
forming regions ( 1 Andersen et al.l 120081 ; IScholz et al.l 120091 ) . 



Bontemps et al.l (120011 ) provide constraints on the stellar mass function in this region, 
complete down to below 0.1 M . From their Fig. 8 we infer a lower limit of ~ 80 stars in 



the g iven mass range. Taking into account an e stimated binar i ty of ~ 29% (IRatzka et al. 



20051 ) would increase this number to about 103. IWilking et al.l ( 120081 ) summarized a list of 
316 L1688 association members with at least a K s -band 2MASS detection, as well as one or 
more signs of association membership, and 169 out of 316 sources have an optical and/or 
IR spectral type estimate. Of these 169, 19 association members are Brown Dwarfs with 
one or more spectral type estimations of M6 or later, while 140 are low mass stars with 
spectral type estimates between G2 and M6 (0.1 < M < 1.0 M Q ). Combining mid- and 
near-infrared as well as X-ray data, these are the most complete estimates for the stellar 
population in p Oph to date. Making the reasonable assumption that the stellar sample is 
more complete than the brown dwarf sample, and adopting a range of 103-140 in number 
of stellar sources, and 19 brown dwarfs, then an upper limit for the ratio low-mass stars to 
BDs is 5-7. Depending on how many BDs the current census is missing, this ratio might be 
as low as ~ 3. As it stands now, the ob served ratio is fully in line with the literature values 
for most other young clusters (3.3-8.5, lAndersen et al.ll2008f). althou gh significantly higher 



than the value we have obtained for NGC1333 (1.5. 



Scholz et al.ll2009l ). Note that the values 



quoted by lAndersen et al.l (120081 ) should be treated as upper limits as well, since they are 
limited to BDs with M > 0.03 M . Thus, there is no reason to believe that p Oph has a 
substellar IMF significantly different from other star forming regions. 
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6. Conclusions 

• Large-scale, optical+near-infrared imaging surveys with Suprime-Cam and MOIRCS 
have been used to create a catalog of sources, reaching completeness limits of 24.2 in 
i-band, 20.6 in J-band, and 17.8 in K s -band. The spatial coverage of our i'+JK s survey 
is 0.171 sq. deg. and covers the L1688 core in i'-, J-, and K s -band. 

• In terms of object masses for members of p Oph, the survey completeness limits cor- 
respond in the i'-band to mass limits of 0.004 - 0.1 M , and in the J-band to mass 
limits of 0.001-0.007 M , for extinction ranges of Ay = 5-15, based on the COND03 
and DUSTY00 evolutionary tracks. 

• From the optical + near-infrared photometry, 309 objects were selected as candidate 
substellar and planetary-mass cluster members. 58 of these objects, and 1 additional 
previously known BD candidate, were targeted for follow-up spectroscopy. Based on 
multi-object spectroscopy in H-band and K s -band, using the water absorption features 
in the H-band, 1 of the 58 new candidates was confirmed as a substellar mass object 
with T cfr = 2500 K. 

• Based on literature and this survey, we identify a sample of 21 spectroscopically con- 
firmed young Brown Dwarf members within our survey area of p Oph. However, this 
sample, based on current existing surveys, is likely highly incomplete, due to the vari- 
able and high extinction in p Oph. 

• From MOIRCS, 2MASS, and Spitzer photometry a sample of 27 sources with mid- 
infrared colour excess and near-infrared colours indicative for substellar mass sources 
with disks are identified. Of these, 11 are previously spectroscopically confirmed brown 
dwarfs, while 16 are newly identified candidates. 

• Based on present day surveys of the stellar and brown dwarf populations, the ratio of 
substellar to stellar sources in p Oph is derived to have an upper limit of 5-7. This is 
in line with other nearby young star forming regions. 
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